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ABSTRACT 

A method is presented to minimize the residuaJ matrix stresses in metal matrix com- 
posites. Fabrication parameters such as temperature and consolidation pressure are opti- 
mized concurrently with the characteristics (ie., modulus , coefficient of thermal expansion, 
strength, and interphase thickness ) of a fiber-matrix interphase. By including the inter- 
phase properties in the optimization of the fabrication process, lower residual stresses are 
achievable. Results for a ultra-high modulus graphite (Pl00)/copper composite show a 
reduction of 21% for the maximum matrix microstress when optimizing the fabrication 
process alone. Concurrent optimization of the fabrication process and interphase proper- 
ties show a 41% decrease in the maximum microstress. Therefore, this optimization method 
demonstrates the capability of reducing residual microstresses by altering the temperature 
and consolidation pressure histories and tailoring the interphase properties for an improved 
composite material. In addition, the results indicate that the consolidation pressures are 
the most important fabrication parameters, and the coefficient of thermal expansion is the 
most critical interphase property. 
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Shear Modulus. 

Inequality constraint. 
Interphase thickness. 

Fiber volume ratio. 

Property. 

Pressure. 

Ultimate Strength. 
Temperature. 

Melting temperature. 

Time. 

Weight coefficient. 
Optimization Parameter. 
Thermal expansion coefficient. 
Mass Density. 

Stress or Microstress. 

Poisson’s ratio. 


Compression. 
Lower bound. 
Tension. 
Upper bound. 


Interphase. 

Fiber. 

Matrix. 

Composite. 

Normal. 

Reference. 

Shear. 

Composite material coordinate system axes. 


2 



1. INTRODUCTION 


Metal-matrix composites (MMC) have generated a high-degree of interest for 
aerospace and structural applications. Despite their high cost, MMCs are potential can- 
didates for specific applications demanding high operational temperatures (from 400 C to 
1100°C), hygrothermal resistance, stability, and peak mechanical performance. A crucial 
problem limiting the use of many MMCs is the high residual (final) thermal microstresses 
developed during the fabrication process, as a result of the large temperature differential 
and the mismatch between the thermal expansion coefficients (CTE) of the fiber and ma- 
trix. The presence of residual microstresses typically degrades the mechanical performance 
of the composite and is primarily responsible for the reported poor thermo-mechanical fa- 
tigue endurance of many MMCs. 

It is desirable, therefore, to explore possible ways to reduce, or alternatively control, 
the development of residual microstresses. One possibility is to use a suitable fiber coating 
as an interphase layer between the fibers and the matrix to reduce the effects of the 
fiber/matrix CTE mismatch [1]. Recent work [2] has also shown that reductions in residual 
stresses may be also attained with optimal combinations of temperature and consolidation 
pressure during fabrication. Hence, it seems advantageous to tailor not only the interphase 
properties, but also the temperature and consolidation pressure profiles of the fabrication 
process concurrently to achieve acceptable residual microstresses. 

This paper presents the development of methodologies for: (1) the simulation of the 
micromechanical nonlinear response of unidirectional MMCs with three material phases 
(fiber, matrix, and interphase) during fabrication; and (2) the concurrent tailoring of the 
interphase layer with optimal fabrication considerations. The objective is to minimize the 
residual matrix microstresses at the end of the fabrication process by optimizing the tem- 
perature and consolidation pressure histories concurrently with the interphase properties 
(modulus, CTE, strength) and other composite parameters (interphase thickness and fiber 
volume ratio (FVR)), while the material integrity throughout the process is ensured. By 
concurrently optimizing the process and the interphase, greater 
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reductions in residual microstresses are achieved when compared to individual opti- 
mization of the interphase and fabrication process. Finally, applications of this method 
are shown for an ultra- high modulus graphite (P100)/copper composite system. 

2. FABRICATION PROCESS 

A typical fabrication cycle, as the one shown in Fig. 1 for PlOO/copper, usually 
involves three phases. During Phase 1 the temperature of the raw materials (rows of fibers 
stacked between ultra-thin foils of matrix, or rows of fibers in metallic powder, or tows 
of graphite fibers preimpregnated in the metallic matrix by the chemical vapor deposition 
method) is elevated near the matrix melting point. Usually, the maximum temperature 
remains lower than the melting temperature of the matrix to reduce chemical reactions at 
the fiber-matrix interphase and matrix oxidation, which requires the application of higher 
consolidation pressure (Phase 2) such that the matrix is diffused and bonded with the fibers 
(superplastic diffusion bonding). Phase 2 is assumed sufficiently long, such that perfect 
bonding is accomplished and pre-existing residual stresses in the matrix are relaxed. In 
the final cool-down phase of the composite (Phase 3) the temperature and consolidation 
pressure are reduced to room conditions. 

Phase 3 is the more critical phase of the fabrication process because significant residual 
microstresses are developed in the matrix and the fibers, as a result of the mismatch be- 
tween the CTE of the fibers and matrix. Mechanical microstresses are also present during 
phase 3 as a result of the consolidation pressure, but they vanish at the end of the fabri- 
cation when the consolidation pressure is removed. Since the build-up of residual stresses 
and the integrity of the composite material are primarily affected by the temperature and 
consolidation pressure histories during Phase 3, it is reasonable to focus the current study 
on this stage of fabrication. 

3. THERMO-MECHANICAL RESPONSE 

The thermo-mechanical response of MMCs during Phase 3 of the fabrication process 
was based on a unified non-linear micromechanics theory developed by Chamis and co- 
workers [3-5]. The theory incorporates, among other factors, three material phases (fiber, 
matrix, and interphase), temperature effects, and mechanical non-linearities. An in-house 
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code (METCAN: METal-matrix Composite ANalyzer) was used to simulate the thermo- 
mechanical response of the MMCs. The basic elements of the method are summarized in 
this section. 

The micromechanics theory assumes the composite microstructure shown in Fig. 2 
consisting of three material phases: the fiber (/), the matrix (m) , and an interphase (d). 
In this paper the interphase represents the fiber coating placed between the fiber and the 
matrix to reduce the build-up of thermal residual stresses in the matrix. Three distinct 
micro-regions are recognized in the composite material [5], as shown in Fig. 2, which are 
identified with subscripts A (matrix), B (matrix-interphase), and C (matrix-interphase- 
fibers). Micro-regions A, B , and C have different average properties, therefore, mechanical 
or thermal loads result in the development of different average stresses. 

The average properties of the composite, as well as, the average properties of rmcrore- 
gions A, 5, and C are provided as explicit functions of the properties and the volume 
ratios (&) of the three constituents. In general, 


Pj — f {P^n ’ Pd* Pfi^fikd) j — A, B,C 


( 1 ) 


where P represents a specific property and subscript j represents either the composite (/) 
or microregions A, B , and C . Subscripts m, d, and / identify the matrix, interphase and 
fiber respectively, and superscript t properties at time t . The theory assumes constant 
average stresses in each microregion, admissible boundary conditions, and equilibrium of 
stresses. The assumption of constant average stress mandates the correction of the fiber 
and interphase volume ratios in the transverse and shear directions. The latter assumption 
has been extensively applied in the development of micromechanics for polymer composites, 
and has provided good accuracy at minimal computational cost. 

In the present paper the nonlinear effects of state variables, such as temperature and 
stress, on the in situ properties of the constituent materials are represented by the following 
form: 


EL 

Poi 


Tvi-T'-rSt-crt 


i = m, d, / 


( 2 ) 
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where: subscript i indicates matrix, interphase, or fibers; subscript o reference conditions; 
subscript M the melting point. Candidate properties for this equation are the moduli 
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(jE 7), Poisson’s ratios^), strengths (5), and CTE (a) of the constituents. The first term 
in the right hand side of eq. 2 represents the temperature effects and the second term 
the mechanical nonlinearity in a stress-strain curve. Time effects have been neglected. All 
exponents in eq. 2 are estimated from correlations with experimental data. 

Because of the non-linear behavior of the material, an incremental procedure was used. 
In this context, the mechanical composite stresses, the temperature, and the microstresses 
at time t -f At are from the respective cumulative quantities at time t and their increments 
during time step At, that is: 

M +A ‘} = «} + W (3-D 

r‘ +ai — T*+T (3.2) 

M (3-3) 

where time superscripts and no superscripts indicate cumulative and incremental quantities 
respectively. The subscript j represents each microregion. 

The incremental microstresses in eq. 3.3 at the different regions of the composite 
induced by the incremental stresses {cr*} and temperature T, are calculated from closed- 
form expressions derived on the previously stated micromechanical assumptions and linear 
elastic response for the homogenized composite and the individual constituent materials 
during time step At. 

= (4) 

The properties of the constituent materials Pj in eq. 4 are calculated from eq. 2, and 
the equivalent composite properties Pf from eqs. 1-2. The highly nonlinear behavior of 
the matrix and the interphase requires either the use of sufficiently small time steps At, 
or combinations of larger time steps and an iterative procedure. The latter technique was 
used in the present paper. 
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4. FABRICATION/INTERPHASE OPTIMIZATION 


The proposed method aims to minimize the residual matrix microstresses by optimiz- 
ing: (1) the consolidation temperature and pressure, and (2) the mechanical characteristics 
of the interphase. Considering the large number of parameters and the complexity of the 
simulation, this may be best accomplished with non-linear mathematical programming 
(NLP). It is recalled that a standard constrained NLP problem involves the minimization 


of an objective function: 


min F( z) 

(5.1) 

subject to constraints in the following form: 

z L < z < z u 

(5.2) 

Q{ z) < 0 

(5.3) 


In the present paper emphasis is placed on the minimization of the matrix stresses 
in region A. In the case of open-die consolidation (ie. application of equal pressure in 
both transverse directions 22 and 33, and no pressure in the longitudinal direction 11), 
only the normal microstresses <T mAn and cr mA2 2 , (where er m A22 = exist in the 

matrix (region A). Among the many possible ways for these stresses to be minimized 
simultaneously, the mini-max formulation, ie. minimize the maximum stress, is proposed 
for its tendency to result in equal minimum stresses. Therefore, the optimal fabrication 
problem is first formulated as the following constrained optimization, 

imn(imx{wio- mAU ,w 2 a mA22 }) ( 6 ) 

subject to upper and lower bounds (5.2) on the optimization vector z. The optimization 
vector includes: (1) the temperatures, consolidation pressures, and times at n p control 
points defining, n p — 1 segments of linear temperature and pressure variations; and (2) 
critical mechanical properties of the candidate interphase in reference conditions, such as 
the modulus, ultimate strength, CTE, and so forth. 

Constraints are also imposed on the matrix, interphase, and fiber microstresses at n, 
time steps in the form of the maximum stress criterion, 

^ ( a mkll )j < )j J ~ ^ n V anc ^ k = A, B (7 .1) 
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i^mk22)j < ( cr mfc22); < (^"mjt22 )j 

j = 1, ..., Tip and k = A, B ,C 

(7.2) 

{^dkll)j < {^ilkll)} < {Sjkll)j 

j = 1 , ...,n p and k = B,C 

(7.3) 

{^dk22)i < ( a dk22)j < (Sjk22)j 

j = 1, ..., rip and k = B,C 

(7.4) 

(Sfcu)j < ( Cr fcn)j < (Sjci i)j 

j = 1 

(7.5) 

{SfC22)j < ( C7 7C22); < {Sj C22 )j 

3 — 1 , n p 

(7.6) 


The superscripts C and T identify compressive and tensile strengths respectively. An 
additional constraint is imposed on the interphase thickness to ensure topological com- 
patibility in the case for square packing of fibers: 




(7.7) 


Whereas, df is the fiber diameter and kf is the FVR. 

The optimization criteria described by eqs, 6-7 are transformed to an equivalent NLP 
compatible formulation (eqs. 5) as follows: 


min(C) 


( 8 . 1 ) 


subject to constraints, 

WlVmAU < C (8.2) 

W20W22 < C (8.3) 

in addition to constraints (5.2) and (7). The objective function £ is an additional design 
variable. The NLP problem described by eqs. 8, 5.2, and 7 is numerically solved with 
the modified feasible directions non-linear programming method [6]. The modified feasible 
directions algorithm performs a direct search within the feasible optimization domain. 
The search direction is estimated from first order sensitivity of the objective function 
and the active constraints. A line search follows along the calculated search direction. 
The implemented algorithm includes an active set strategy, ie., only the constraints near 
violation are included in the search, thus allowing the efficient handling of the large number 
of constraints defined by eqs. (7). 
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5. APPLICATION AND DISCUSSION 


An ultra-high modulus graphite (P100)/copper unidirectional MMC was used to test 
the optimization method for the fabrication process and the concurrent optimization of 
the process and interphase properties. Representative constituent properties of the com- 
posite system at reference conditions (21° C, 0 Pa) are shown in Table 1. The Materials 
Division of NASA Lewis Research Center provided the current fabrication process for the 
PlOO/copper. Only the cool-down phase of the fabrication process will be simulated dur- 
ing the optimization and thermo-mechanical response of the MMC. The development of 
residual stresses and the integrity of the composite material are primarily affected by the 
temperature and pressure histories at this phase. Initial interphase properties were as- 
sumed equivalent to the matrix properties. In addition to the material properties, the 
initial interphase thickness was 12% of the fiber diameter and the FVR of the composite 
system was 40%. 

The cool-down phase was subdivided into four increments of linearly varying tem- 
perature and pressure. Stress constraints were imposed at five evenly spaced time in- 
tervals in each linear segment. In this manner, twenty constraints were introduced for 
each microstress inequality described in eq. 7 , and where the interphase was optimized 
one additional constraint on the thickness will be added. The weighing coefficients were: 
w\ = W 2 = 1. The temperatures, pressures, and times at the starting and final points of 
the four linear segments were used as optimization parameters. The temperature at the 
beginning of the cool-down phase was held constant and equal to the respective temper- 
ature of the current processes, and the final pressure was set equal to zero. Shown below 
are the upper and lower bounds imposed on the optimization method in accordance with 
eq. 5.2. 

Fabrication Process Parameters: 


T 0 <T< T m 

(9.1) 

0 < p < 345 MPa 

(9.2) 

lOsec < i < 18000sec 

(9.3) 
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Interphase Property: 


34 .6GPa <E d < 220.8 GPa 

(9.4) 

1.69 cm/ cm/°C < a d < 67.8 cm/cm/°C 

(9.5) 

ZA.hMPa <S d < 414.0 MPa 

(9.6) 

Micromechanical Parameters: 


0.05 <k d < 0.15 

(9.7) 

0.05 < k f < 0.55 

(9.8) 


Figure 3 shows the current and the optimum fabrication processes, Case 1 (the fab- 
rication process optimized alone) and Case 2 (concurrent optimization of the fabrication 
process and interphase properties), for the PlOO/copper MMC. Both optimized processes, 
Cases 1 and 2, follow similar patterns during the cool-down phase. Compared to the cur- 
rent process, two significant differences exist that lead to the reduction of the final residual 
matrix microstresses. First the optimized temperature histories in Fig. 3a decrease more 
rapidly to room temperature and are held constant until the end of the processes; secondly 
as shown in Fig. 3b, the predicted optimal consolidation pressures gradually increase as 
the consolidation temperature drops, reaching significantly higher values than the pressure 
of the current process and finally dropping to zero. 

More interestingly, the temperature drop takes place when the pressure is high, indi- 
cating that the thermal stresses are forced to develop when the matrix and interphase are 
highly nonlinear and near “plastic”, hence high strains result in low stresses. The pressure 
is removed when the temperature reaches room values as it does not contribute further. 
This illustrates the importance of the consolidation pressure. Apparently the values of the 
consolidation pressure were limited by the current strength of the constituents, as indicated 
by the observed active constraints in eq. 7. Shown in Fig. 4 are the predicted buildups 
of microstresses 411 and a m A 22- The final (residual) microstress <r m Ai 1 Case 1 de- 
creased by 21% compared to the respective microstress value of the current process, in 
comparison, Case 2 had a 41% reduction for CmA 11 • The additional reductions in Case 
2 were attributed to the interphase optimization. Microstress <TmA22 °f the Case 1 was 
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nearly equivalent to the current process, however, the final microstress 22 decreased 
by 24% because of optimized interphase properties. 

The optimization of the interphase properties between the fiber and matrix, refer to 
Fig. 2, resulted in an optimum interphase with compatible properties (Table 2) to the 
fabrication process and fiber. The fact that the initial interphase material was matrix 
indicates the need for an interphase between the fibers/matrix. As seen in Table 2, all 
interphase properties increased during the optimization. Specifically, the modulus and 
strength increased slightly, while the CTE had a larger increase. As indicated by the 
results, the CTE was the most critical interphase property. The interphase thickness and 
FVR also increased. 

6. SUMMARY 

A method was proposed for optimizing the fabrication process and interphase prop- 
erties concurrently for unidirectional metal-matrix composites. The response of the fabri- 
cated MMC was simulated based on nonlinear micromechanics and the NLP problem was 
numerically solved with the modified feasible directions nonlinear programming method. 
Reduction in residual microstresses for the fabricated composite were optimized by varying 
temperature and pressure histories along with the interphase properties during the cool- 
down phase. An in-house research code has been developed incorporating this method. 

Case studies were performed on ultra-high modulus graphite (P100)/copper compos- 
ite. The predicted optimum process for the fabrication process reduced the maximum final 
microstress hy 21%, by favorably optimizing the nonlinear in situ matrix behavior. By in- 
cluding the interphase properties in the optimization method a further reduction, 41%, 
in the final maximum microstress was achieved illustrating the importance of the inter- 
phase optimization. Finally, the results indicate that consolidation pressures are the most 
important fabrication parameters, and the CTE is the most critical interphase property. 
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(a) Temperature (b) Consolidation Pressure 


Figure 1. - Typical Processing Cycle for Graphite/Copper Composite 



A: Matrix 

B: Matrix and interphase 
C: Matrix, Fiber, and Interphase 
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Fiber 


Figure 2. - Material Microregions in a Representative MMC Cell 
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(a) 


(b) 


Figure 3. - Optimum and Current Cool-Down Phases for PlOO/Copper 


Current 



(a) Longitudinal Stress (b) Transverse Stress 

Figure 4. - Matrix Microstresses Developed During the Cool-Down Phase of PlOO/Copper 
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